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1. INTRODUCTION 

The grasshopper algorithm (GOA) is a new optimization algorithm developed in 2017 [1] to simulate 
the behavior of grasshoppers in nature. It has been used in different applications. GOA is formulated to 
simulate the behavior of grasshoppers. Different variants of GOA have been suggested. Binary GOA 
schemes are developed [2], [3]. In several works authors use chaotic function to enhanced GOA operation for 
different applications [4]-[12]. Results illustrated the good quality of the developed algorithms compared to 
GOA and other algorithms. An enhanced GOA scheme using a combination of levy flight with GOA 
(LGOA) is developed in [13]. The performance of LGOA was tested and compared with PSO and the basic 
GOA. Results illustrated the good quality of the developed algorithm compared to the basic GOA and PSO 
algorithm. Authors in [14] presented a dynamic GOA using dynamic weight and random jumping concept. It 
was illustrated that the suggested method has good properties. An adaptive GOA is developed to enhance the 
performance of GOA and used for optimization in different applications [15]-[17]. Results show good 
achievement for the proposed method. GOA is used in [18] for tuning the parameters of the Support Vector 
Machine. A dynamic GOA is developed and used for the unmanned aerial vehicles in [19]. The 
multiobjective GOA is used for solving the benchmark problems in [20]. 

An improved grasshopper using levy flight is presented in [21], [22]. Ewees et al. [23] an 
opposition-based learning approach to the GOA is presented. A fuzzy logic is used with GOA for electric 
vehicle charging sites in [24]. The new developed fuzzy GOA was tested on distribution networks and results 
show good performance. A GOA with fuzzy control scheme to track the maximum power point is proposed 
in [25]. GOA is utilized for tuning membership functions of fuzzy logic controller. It was illustrated that 
good performance was obtained. 
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In general the meta-heuristic algorithms use a random initial population for the considered problem, 
find solutions using the objective function. Algorithms enhance the solutions iteratively until satisfying the 
stopping condition and find the global optimum in different applications [26]-[29]. These algorithms have the 
following advantages: they are simple and easy to apply. Independent of derivation. All algorithms have one 
common feature in the searching process for the optimal solution: exploitation and exploration. The 
algorithm in the exploration, go through the global search space. While in the exploitation local search is 
performed. 

The GOA work well and can find good solutions in a reasonable time. In general the GOA have 
linear comfort zone parameter which is weak in balancing exploitation and exploration. Also it may fall into 
local optimum. In [30], [31] authors presented different attempts for the improvement of GOA algorithm. 

This paper proposed a new modified GOA. The major contributions of this work is to replace the 
linear comfort zone parameter by a nonlinear function which results in balancing exploitation and exploration 
and improve the achievement. This article is arranged as: The original GOA and the new MGOA are 
presented in sections 2 and 3 respectively. Results are shown in section 4. Section 5 contains the conclusions. 


2. GRASSHOPPER ALGORITHM 

The GOA is a new metaheuristic algorithm based on actions of long and sudden transitions of 
grasshoppers in groups. Searching for food sources is an important action of grasshoppers [1]. The searching 
process has two stages: exploration and exploitation. The far-reaching and sudden transitions of grasshoppers 
simulate the exploration, and the local transitions in searching for food locations simulate the second stage. 

Pushing and pulling forces between the grasshoppers are described by a mathematical model. 
Pushing forces allow grasshoppers to go in the search area, whereas pulling forces help them in the 
exploitation of the expected areas. The convergence factor in GOA algorithm decrease the comfort zone to 
keep exploration and exploitation in balance. Finally, the current best solution was treated as a target to be 
improved. The mathematical model is as [1]: 


xj = S; +G+Ai (1) 
where S; is the social relation: 

Si = Xita jai S(dij) dij (2) 
where, 


N: the total number of grasshoppers 
dij: the space between grasshopper i and grasshopper j 


dij = |x; — x;| (3) 
u (4) 


xi: is the location of the grasshopper I, G in (1) is the force of gravity, defined as: 
G=gé (5) 


where g is a constant and eis a unit vector in the direction of the center of earth and A is the direction of 
wind, defined as: 


A= uey (6) 


u is a fixed factor and e_w is a unit vector in the wind direction 
The model described in (1) could be redefined as: 


x; = DperjaiS (dij) dij —gétuey (7) 
where, 


S(r) = fe" se (8) 
where, 
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f: attraction constant 
l: attractive length scale 

The effects of wind and gravity are small when compared with the effects of relationships between 
grasshoppers, the model can be redefined as (9): 


hb-sb 


xi = c Œa jet E s( dij) dij) + Ad (9) 


where hb is the higher value and sb is the smaller value in the considered space, Ad is the current best 
solution, c is a comfort zone parameter which is described by: 


c=- iter F (10) 
where c; is the larger value (c=1), cs is the smaller value (c;=0.0004), iter means the running iteration, and 
Maxiter is the highest number of iterations. 

GOA is equipped with a convergence parameter that iteratively decreases the comfort zone to 
reconcile the exploration and exploitation and skip from local optima. It is an important coefficient and can 
be considered as a convergence parameter. The linear comfort zone parameter proposed in (9) may be not 
able to achieve the balancing between the exploration and exploitation. In order to enhance the achievement 
of the GOA a nonlinear comfort zone parameter is suggested to replace the linear parameter. The comfort 
zone parameter is suggested here to be as (c™, m=2, 3,...). In the case of linear parameter, the value is 
linearly decreasing with iterations while in the modified algorithm here the value of the parameter is 
proportional to the square of the number of iterations (for m=2). Figure 1 shows values of the comfort zone 
parameter against iterations for different values of m. The value (m=1) means the parameter in the original 
algorithm, As can be seen the effect of changing m on the curve shape. They are giving different searching 
mechanisms. 


convergene factor 


iteration 


Figure 1. Comfort zone parameter for different values of m 


3. MODIFIED GRASSHOPPER OPTIMIZATION ALGORITHM 
In GOA a linear convergence parameter is used. The linear comfort zone parameter is not able to 
achieve the balancing between the exploration and exploitation In order to enhance the performance of the 
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GOA a nonlinear parameter is suggested to replace the linear parameter. In the modified algorithm the 
comfort zone parameter is chosen as: 


Maxiter 


x m 
c= |c - iter -2] (11) 
where m=2, 3,... 


Clearly the comfort zone parameter c directly affects the change amount in (9). If the parameter 
value is large the change amount is also large. The magnitude of change must be big in the exploration to 
keep going through the global optimum and skip from local optimum. The amount of change must be small 
in the exploitation to make the algorithm smoothly go in the global optimum. In original GOA, magnitude of 
the comfort zone parameter reduces linearly in the two stages. Also, the amount of change in (9) is fixed and 
cannot be arranged according to the requirements of the two stages, which reduces the efficiency of the 
algorithm. In the new modified algorithm the step size varies with the changing of iteration. Also the step 
size depends on the value of m. Each value of m give different step size and different pattern of variation as 
shown in Figure 1. The step size is large at the beginning points which enhance the exploration phase to 
ensure exploring the global optimum and avoid going to local optimum. Finally the step size decreases which 
increase the efficiency of exploitation by converging to the global optimum gradually. In (9) the outer c is 
used to reduce the transitions near the target. The inner c reduces the three zones. The initial population is 
generated using random numbers by the following equation: 


X i = Xmin t rand[0,1] (Xmax =< Xmin) (12) 
where Xmin and Xmax are smaller and larger values of parameter x ; respectively. 


Algorithm (The modified grasshopper optimization algorithm) 


: Select the values c), cs 

: Select the size of population N and maximum number of iterations Maxiter 
: Find Random initial population xi(i = 1,..N) using (12) 

: Determine the fitness f (xi) 

: Let H=Best result 

: While stop is not satisfied do 

: Select the type of the nonlinear convergence parameter by setting the value of m 
: Update the convergence parameter c using equation (11) 

: Fori=1 to N do 

10: Find the normalized distance in the range [1,4] 

11:Use equation (12) to find position of current grasshopper 

12: Keep the current grasshopper inside the boundaries 

13: EndFor 

14: Change H when the solution is better 

15: iter=iter+1 

16: Endwhile 

17: keep the best solution as H 


COMAAINDMNBHWN HE 


4. RESULTS AND DISCUSSIONS 

In order to test the proposed modification. The benchmark used in the GOA original algorithm is 
adopted as a test environment. It consists of different kinds of functions such as unimodal and multimodal. 
The mathematical description and properties of these functions are illustrated in Table 1 and Table 2. Results 
of the modified GOA method have been compared to results of original GOA. All programs were written in 
MATLAB on a PC with the windows 10 and the RAM 8GB. In all test functions for the unimodal and 
multimodal 30 search agents over 500 iterations are used Population size is 100 for each function. In the 
experiment 10 attempts were performed using random initial values to find all results, the Mean fitness. 
standard deviation (SD), the best fitness and worst fitness. Results for the unimodal functions F1-F7 are 
presented in Table 3 (see in appendix). One can notice that the new MGOA show an excellent results in F1, 
F2, F4, F6 and F7 functions for m=2, 3 and 4. Results for m=4 are the best and indicate that the exploitation 
is improved in the modified algorithms. The main reason for the improved performance of MGOA is the 
replacement of the linear comfort zone by the proposed nonlinear scheme. Table 4 (see in appendix) lists the 
results of multimodal test functions F8-F13. The new MGOA works better in almost all test functions. 
Results illustrate that the new MGOA has good exploration properties. 
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Table 1. Unimodal functions 


Function Dim Range Ími 
n 
F1(x) = > x? 30 [-100-100 0 
i=1 
n n 
F2(x) = x+] | Fal 30 [-100-100 0 
m i=1 phe 
F3(x) = 2 > x; ) 30 [-100-100 0 
=1 yt 
F4(x) = Max, flx; | 1<i<n 30 [-100-100 0 
F5(x) = y [100(xi1 — x7)? + (x; — 1)?] 30 [-100-100 0 
i=1 
n 
F6(x) = 3 (lx; + 0.5])? 30 [-100-100 0 
i=1 
n 
F7(x) = > ix} + random[0,1) 30 [-100-100 0 
i=1 


Table 2. Multimodal functions 


Function Dim Range fmin 
p a E * 
F8(x) = 3 -xsin D 30 [-500-500] 316*d 
i=1 
F9(x) = [x? —10cos(27x;) + 10 30 [-5.12- 5.12] 0 


a 30 [-32-32] 0 
F10(x) = —20 exp| —0.2 any =) 


— exp (A) cose 1x;)) +20+e 


1 j a f 30 [-600-600] 0 
2 (2S 2 ; 
F11(x) = G00) Dace xi | [cs (x;/Vi)+1 
T. ; mad 30 -50-50 0 
F12(x) = @A0sin@ny,) ++) O1 = DPA + 10 sin? Ya) Daas 
tOn- DOn -1+ u(x, 10.100.4) 
A i=1 
+ 10.100. 4 
> u0 ) 
Yi= 1+ (x: + 3) 
k(x,;-a)™x>a 
u(x;, a. k,m) -| o-a<x>a 
k(—x; — a)" x < —a 
F13(x) = 0. 1fsin? (37x) +] 30 [-50-50] 0 


+ 2 a — 1)? [1 + sin? (3x + 1)] 
+(x, - T + sin? (31x, + 1)]} 
$ X uas 5,100, 4) 


To show that results were not obtained by chance, the Wilcoxon ranksum test is used. If the 
probability (p) is lower than 0.05 then the difference between the two compared algorithms is significant. For 
this statistical test, all versions of the proposed algorithm in each test function is compared with the original 
algorithm. From Table 5 (see in appendix) it can be seen that the p-values between the original GOA and the 
proposed modified compared algorithms for all values of m in all test functions are less than 0.05. Clearly 
results illustrate the significant improvements in our modified versions of GOA. 

In general the proposed modification results in good improvement for all values of m=2, 3 and 4 
when compared with that of m=1 (original GOA algorithm). For more clarity, the best solutions were marked 
in bold. It can be seen that the performance of MGOA is better on almost all functions. 


5. CONCLUSIONS 

In this work, a novel version of GOA name MGOA was developed. A nonlinear comfort zone 
parameter is proposed to replace the linear parameter in the first GOA. The MGOA is evaluated using the 
functions presented by GOA researchers and compared with the original GOA. Results show that the new 
MGOA has an excellent exploitation and exploration capabilities and show that the new modified MGOA is 
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effective and able to determine the optimal solutions in all functions. Experimental results indicate that the 
proposed MGOA outperforms original GOA, presenting a higher accuracy, faster convergence, and stronger 
stability. The future work will focus on two ideas. In the first different initial population generation strategies 
would be examined and in the second idea, further studies on the selection of the comfort zone parameter 
would be performed. Also we would find additional methods for enhancing the performance of GOA and use 
it in engineering applications. 


APPENDIX 
Table 3. Results of unimodal functions 
Function M Mean SD Best Worst 
F1 1 2.1678e-06 6.6460e-06 3.4741e-07 6.1939e-06 
2 1.7792e-12 4.8005e-12 2.8956e-13 4.517e-12 
3 8.0000e-16 2.2295e-15 1.1104e-17 1.8334e-15 
4 5.0648e-21 2.0068e-20 2.2751e-22 1.6396e-20 
F2 1 0.0242 0.1153 0.00032235 0.097916 
2 8.1921e-07 7.1790e-07 5.4949e-07 1.0825e-06 
3 2.1116e-08 5.9596e-08 1.707e-09 4.7223e-08 
4 4.1621e-10 2.3132e-09 1.8863e-12 1.8777e-09 
F3 1 0.0589 0.0415 0.030205 0.0787536 
2 0.0082 0.0318 9.8936e-15 0.021164 
3 6.791 1e-05 4.2910e-04 5.9138e-10 2.0549e-5 
4 5.4672e-07 5.1536e-06 5.5631 e-12 5.4358e-6 
F4 1 0.0156 0.0632 0.0009637 0.054691 
2 2.2260e-06 6.0913e-06 3.1023e-7 4.5472e-6 
3 6.791 1e-05 4.2910e-04 1.8964e-8 0.0003393 
4 1.3001e-07 1 8.0817e-07 6.0546e-10 6.4113e-7 
F5 1 207.4406 778.2415 5.285 578.5194 
2 6.9545 4.5632 5.285 8.9299 
3 34.8973 105.5537 6.3543 71.2648 
4 7.5222 5.1641 4.1546 9.6059 
F6 1 1.1800e-06 2.3602e-068 4.8463e-7 2.4919e-6 
2 2.9441e-12 4.7629e-12 1.8766e-12 5.2245e-12 
3 2.8883e-17 3.2975e-16 1.7654e-17 5.1127e-17 
4 3.7280e-21 6.5789e-21 2.7154e-21 5.9828e-21 
F7 1 1.8255e-18 6.7801e-18 4.2132e-19 6.0694e-18 
2 2.8716e-30 2.7632e-30 1.128e-30 3.3738e-30 
3 2.6432e-39 1.0193e-38 3.3218e-44 6.5908e-39 
4 5.2546e-48 3.2254e-47 4.4555e-50 2.5652e-47 
Table 4. Results of multimodal functions 
Function M Mean SD best Worst 
F8 1 -6.9539e+03 1.6211e+03 -7646.987 -6327.604 
2 -8.2242e+03 1.1811e+03 -8640.5274 -7693.292 
3 -8.1603e+03 412.8873 -8302.5009 -8004.5934 
4 -6.1812e+03 1.4696e+04 -7819.1874 -7592.2661 
F9 1 26.1436 32.0982 9.9496 41.7881 
2 19.8991 27.3206 10.9445 39.7983 
3 24.9734 37.7938 4.9748 51.7376 
4 16.4168 10.8309 8.9546 19.8991 
F10 1 4.5754e-04 4.6372e-04 0.00026956 0.0006719 
2 6.5230e-07 6.0253e-07 3.6482e-7 9.4439e-7 
3 5.0896e-09 1.6358e-08 1.6549e-9 1.5182e-8 
4 1.4734e-11 1.4193e-11 3.5039e-12 4.1926e-11 
F11 1 0.2551 0.1876 0.15587 0.31061 
2 0.2254 0.0644 0.19919 0.24359 
3 0.1783 0.2046 0.126105 0.30469 
4 0.1595 0.1737 0.076279 0.2044 
F12 1 1.3859e-04 5.2703e-04 3.7083e-7 0.0003427 
2 1.0131e-06 4.1388e-06 3.459e-10 3.459e-6 
3 2.4014e-09 6.4920e-09 9.842e-12 4.6537e-9 
4 1.5215e-11 5.4786e-11 1.7518e-15 4. 8984e-11 
F13 1 0.0044 0.0171 6.453 e-6 0.011068 
2 4.2959e-08 1.4298e-07 5.7973e-9 9.8332e-8 
3 1.0992e-09 3.8803e-09 1. 245e-10 3.1829e-9 
4 3.2598e-12 1.2862e-11 1.0399e-13 1.0122e-11 
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Table 5. Values of probability from Wilcoxon ranksum test 


Function/m 2 3 4 
Fl 1.6494e-04 1.7265e-04 1.6494e-04 
F2 1.6118e-04 1.6118e-04 1.7364e-04 
F3 1.5748e-04 1.5748e-04 1.1203e-04 
F4 1.5748e-04 1.6494e-04 1.6494e-04 
F5 0.0107 0.1023 0.0206 
F6 1.6494e-04 1.6494e-04 1.6494e-04 
F7 1.4332e-04 1.4332e-04 1.5295e-04 
F8 1.4851e-04 1.4851e-04 0.00196 
F9 0.1612 0.8205 0.0409 
F10 1.6780e-04 1.6780e-04 1.6494e-04 
F11 0.0738 0.0108 0.0071 
F12 0.0043 1.5748e-04 1.5748e-04 
F13 1.5748e-04 1.5748e-04 1.5748e-04 
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